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film thickness; H = h/C 
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coordinates located in centre plane of bearing at intersection of the line of 
centres O b O j with the bearing surface as shown in Figure 1; Z - Z / L 

X coordinate of the j'th oil hole 
series coefficients defined by eqn. (20) 

angular location of the point A from the line of centres along the bearing 
surface as shown in Figure 1; 0 = X/R 
angular location of the j'th oil hole; 0j * Xj /R 
mean viscosity of lubricant 

polar coordinates with origin at centre of j’th oil hole as shown in Figure 3 
unbalance eccentricity in Figure 1 

time dependent multiplier of Pj 
defined by eqn. (19) 

angular displacement of the line of centres from the vertical direction as shown 
in Figure 1. 

a constant; T] = 1 for journal bearings, T] — 0 for dampers 

W b +(Oj 

angular velocity of bearing and journal respectively 
L/dj 

denotes differentiation with respect to time x 

INTRODUCTION 

Because of their compact size and wide range of damping capabiity, squeeze film dampers 
are frequently used for stabilisation and vibration attenuation of rotating machinery. Oil 
feed is generally via a central circumferential groove, effectively splitting the damper width 
L into two dampers, each with axial land length of L/2. However, such a design effectively 
reduces the damping capability of the damper by a factor of four. To recover this lost 
damping capability, one can dispense entirely with the circumferential groove feed and 
introduce all the lubricant from either of the bearing ends. Such a concept requires an axial 
seal to ensure that the lubricant does not bypass the damper clearance space. A practically 
more reliable alternative is to introduce the lubricant via one or more (usually two or four) 
oil holes in the central axial plane of the damper, i.e. replacing the circumferential groove 
w ith a number of oil holes. While this arrangement ensures that al l the lubricant passes 
through the damper clearance, the unbalance response analysis for even the simplest case of 
a squeeze film damped centrally preloaded horizontal rigid rotor or a squeeze film damped 
vertical rotor has become significantly more complicated. Whereas with the end feed 
dampers, one can assume circular orbit response in the steady state, such is no longer the 
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■ w,,h 0,1 hole feed dampers where the steady state response generally needs to be 
obtained by a transient solution. Transient analyses rely heavily „„ approximations such as 

bear '" 8 Wr0xima,i °” or compilation of appropriate tabular data, to be 

anD> host Uy H eaS,b ' e ’ ' he Ia,ler aPPr “ Ch be '" S War,an,ed 81 relalivel !' h iSh aspect ratios 
(17 0.-5) and eccentricity ratios (c > 0.7) (Hahn, 1989). Such a tabular approach is 

xpeeted to be necessary with oil hole feed, so the major emphasis of this paper centres on 

ho* to ob, am satisfactory tabular date and thereby overcome the added compLon due to 
the introduction of oil hole feed. 



figure 1 : Simple circular bore hydrodynamic bearing or squeeze film damper. 

SOLUTION STRATEGY 

Figure 1 depicts a simple circular bore hydrodynamic journal bearing or squeeze film 
damper tf the journal does no, rotate (i.e. if co, - 0). Figure 2 shows the journal oil feed 

pressure boundary conditions which now, in addition to spec, fled inlet and outlet pressures 
Pi and Po , allow for „ oil feed holes a, the centre plane of the bearing (i.e, along Z = 0) 
Fbr the usual assumptions of laminar flow, Newtonian fluid, cons, an, fluid properties 

* ‘g. ' ner " a f ° rCe - b0dy f0rce and c-'vature effecls, no slip a, bearing surfaces no 
ax 1 3 Va ' ,al, ° n ” f " m ,hicknKS ' ““*»>' pressure across the film and validity of ordei of 
ntegmtude cons, derations, one can derive , he Reynolds equation as (Pinkus and Sternlieht, 
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Figure 2: Bearing or damper with end feed and end drain pressure boundary conditions. 

The boundary conditions pertaining to eqn. (1), as shown in Figure 2, are: 

p = Pi at Z = - ^ (3) 

P - Po at Z “ Yl ■ 

These boundary conditions describe end feed conditions. Additionally, there are now n 
circumferentially equispaced oil holes in the mid circumferential plane at Z = 0, supplied 
from an oil chamber at supply pressure p s . Finally, the pressure distribution is periodic in 
the circumferential (i.e. 0) direction with period 2ji. 


Then, in non-dimensional form 
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where H = (1 + t cos6) . 

Ignoring for the time being the problem of cavitation, the solution of eqn. (3) may be 
written as 
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respectively, with boundary conditions P s - P c - 0 at Z - ± 
and P s (6) = P s (0 + 2k); P c (e) - P c ( 0 + 2k); 

and wherePj(j - 1,2.... n) arc given by the solutions of: 



L 2 d_ 
R 2 dQ 
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Pj - 1 at Z - 0,0 = 0j. 


(10) 

( 11 ) 


( 12 ) 

(13) 

(14) 


In eqns. (6) and (7), the Aj's are time dependent multipliers so that at any instant of time, 
at Z - 0 and 0 = 0J 

A lPl + A 2 P 2 +A n P n -P-P*-P oj _p* ( 15) 

Strictly speaking, this formulation implies that the pressure will be P 0 j at the centre of the 
oil hole but may be somewhat different from this elsewhere; whereas in fact the pressure 
should equal P oj everywhere within the oil hole region. The error due to this assumption 

will increase as the hole diameter dj increases but owing to the near parabolic axial profiles 
of P s and P c , is expected to have insignificant effect for X s 6. 


If there is flow restriction between the supply pressure P ss and the pressure at the hole at 
0-0j, c.g. a capillary or orifice, one needs in general to balance the flow through the 
restrictor to the flow out of the oil hole. To cater for such an eventuality, an iteration is 
introduced whereby the hole pressures P 0 j are guessed and subsequently modified until 

there is flow balance. Eqn. set (15) then constitutes a set of n simultaneous linear 
equations in the n unknowns * h ...A n . If the restrictors offer sufficiently small 
resistance to the flows, the pressure drop across them will be negligible, whereupon all 
P 0 j - P^ and the need for iteration is avoided. In any case, eqn. (6) gives, for this instant 

of time, the full pressure distribution around the damper. Such an approach allows for fluid 
to be pumped out as well as sucked in at the holes. Such pump out would in practice be 
avoided by using a check valve. In such a case, whereas the pressure can never fall below 
P ss at any oil hole, any respective hole can be ignored once P* exceeds P ss , thereby 
reducing the number of 'effective' oil holes n. 
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The major computational problem involves obtaining the pressure distributions 
Pc Pn and P; at each instant of time. It is computationally prohibitive to solve 

numerically the partial differential equations (8), (9) and (12) at each instant of time, so a 
data tabulation/interpolation strategy is adopted. As may be seen from eqns. 
(8) - (11), P s and P c are, for a given aspect ratio L/R, functions of 0,Z and E . Three 
dimensional arrays for P s and P c with sufficient elements to allow accurate interpolation are 
both impractical and unnecessary for, as shown in Rezvani and Hahn (1992), by assuming a 
parabolic shape for the Pg and P^ in the axial direction, two dimensional arrays which are 
functions of 6 and e only will suffice with adequate accuracy for aspect ratios as high as 
L/D = 1, and eccentricity ratios as high as e = 0.9, or even higher. 


The same cannot, unfortunately, be said for any of the Pj pressure distributions, since they 


are a function of the j'th oil hole location 0j with respect to the line of centres, which 



Figure 3: Laid out journal bearing or damper with oil hole at 0 — 0j. 


location changes with time. Figure 3 depicts the laid out journal bearing or damper with the 
instantaneous location of the oil hole at 0j. Attention can be restricted to the solution for Pj 

for a single oil hole, since the value of 0 j is arbitrary and can therefore apply to any of the 
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oil holes. Thus a sufficient number of Pj distributions need to be stored, to enable Pj to be 

found with sufficient accuracy at any time by interpolation. Additionally, the Pj are also 

functions of 0, Z and e suggesting the need for four dimensional arrays. Nor can they be 
found by approximate or numerical solution as easily as the P s and P c , owing to the finite 
size of the oil holes, and indeed the size of the axial width of the bearing land relative to the 
hole diameter \j is a parameter which is likely to have a significant effect. For finite sized 

oil holes, numerical solution of eqn. (12) using finite differences is awkward, owing to the 
mix of curved and straight boundaries. Thus, Marmol and Vance (1978) assume the oil 
holes to be smaller than the finite difference mesh. Not that a finite difference solution 
would solve the tabulation problem since at any instant of time, the Pj's would be functions 

not only of 0,Z and e but also of 0j . Nor has a satisfactory approximate solution approach 

been found. Rather, a cosine power scries solution approach, described below, has been 
found to be satisfactory and promises to provide manageable tabulation data. 

Once all the pressure distribution components have been found, the full pressure distribution 
P can be found from eqn. (6). Cavitation can then be approximated by setting P = 0 when 
P < 0, a commonly adopted assumption in transient analysis, as the insistence of zero 
pressure gradient at the cavitation boundary would unduly prolong the time required to 
evaluate the pressure distribution, and hence, the fluid film forces, whereas the resulting 
gain in accuracy is expected to be minimal. 


o INTERIOR POINT 



Figure 4: Solution domain for Pj (M = 10, N - 16). 
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COSINE POWER SERIES SOLUTION FOR Pj 
Required is the solution of 
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(16) 


over the domain in Figure 4, where Pj — 0 at Z « ±L/2,Pj — 1 around the circumference 
of a circle of radius aj centred at Z - 0,0 - 0j. 


The circular boundary formed by the oil hole, suggests the reformulation of eqn. (16) in 
terms of polar coordinates (r, xp). With origin of these coordinates at the centre of the oil 
hole, eqn. (16) may be written as: 

d 2 Pj 10P^ ! d 2p j | j dPj dlnH | 1 ap j 

dr 2 + r dr r 2 chp 2 v 5r r 2 dxp , 

where H = 1 + ecos0 - 1 + ecosxj 

and where Xj “ ®j +-^-C°stp = 0j + ^-^cosip (19) 

Thus, Xj is a function of r,xp,Xj,L/D and 0j. 


(17) 

(18) 


Hence, for a given bearing with given L/D and Xj, the Pj at any instant of time will depend 
on eand0j. Generally, all n holes would have equal diameters, so that henceforth 
Xj-X. 

In view of the symmetry about Z = 0, a cosine power series solution approach is adopted. 
Letting 

P j"S 2«ij r 1 c °s(j - l)xp (20) 

i=l j=l 

where the integers M and N are chosen sufficiently large to assure convergence. 
Substitution into eqn. (17) and evaluation of eqn. (17) at selected points within the solution 
domain and along the curved, the 'side' and the 'end' boundaries results in MN simultaneous 
equations, for the MN unknowns ajj. 
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RESULTS AND DISCUSSION 


The potnts selected were along rays as shown in Fig. 4. Each ray was subdivided into 9 
equal segments between the oil hole and an outer boundary, so that M = 10 in Eqn 00) In 
thts regard, since the precise boundaries to the left and right are unspecified and may be 
assumed to extend indefinitely, the solution is restricted to finite distances ±AXj on either 

Side of Xj with the value of Pj set to 0 a, these boundaries. To ensure that the solution 

was virtually asymptotic a. eAXj.AXj is increased until no readily discernible difference 
occurred in the respective solutions. 


Using 16 such equispaced ’rays', so that N = 16 in eqn. (20), resulted in 160 linear 
stmul, ancons equations fo, the 160 values of a,,. This was give sufficient accuracy fo, 

the parameter combinations tested. Also, definition of the side boundaries by selecting a 
AXj of four hole diameters was generally found to be adequate for X = 6. Since eqn. (16) 

is simply Laplace’s equation if H is constant, the accuracy of this solution procedure (for 
constant H) was able to be tested against an existing single series solution procedure 
desertbed tn the literature (Takeuli and Sekiya, 1968), and was found to be in agreement. 


To illustrate the variation of the Pj with t and 6j. Figs. 5 to 10 give the solutions for s = 0.5 

and 0.9, with fl, - 0°,90° and 180° for a bearing with IVD = 0.5 and X = 6. As expected, 

flj 0 and 180 , (Figs. 5, 7, 8 and 10) the solution is symmetrical about the vertical 
centreline, but the pressure contours are skewed for the other 0, Maximum skewness 

would occur around 6j - 90°, and for s = 0.9 (Fig. 9) one would probably need to further 
increase M, N and the AXj to improve the accuracy of the small values of Pj. 


Since 6, is solution dependent, data for a sufficient number of discrete 8j values will need 
to be provided to allow for interpolation in a transient programme. A comparison of Figs. 5 
and 8, 6 and 9, and 7 and 10 illustrates the variability of the solutions with e. Again since e 
ts solution dependent, data for a sufficient number of discrete e values will need to be 
provided to allow for interpolation. However, rather than store the Pj pressure distributions 

fora, angeofvalues of e and0j.it will be sufficient tostore the ctj coefficients, thereby 
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Figure 10: Solution for Pj for X. - 6, E - O.9,0j - 180 . 

pressure contributions from other oil holes is negligible until the number of oil holes is so 
large that neighbouring centres are within 5 oil hole diameters of each other, (e.g. for X. = 6 
and L/D = 0.5, such interaction could be neglected provided there were fewer then about 8 
oil holes). Hence, in most practical situations, it is expected that such interaction is 
negligible, so that eqn. (15) simplifies to: 



AjPj - p - p ‘ 


r ss 


( 21 ) 

at Z-OandB-Bj. This ability to treat oil holes independently also considerably 

simplifies the iteration introduced in the case of a flow restriction between supply reservoir 
and the hole exit. 


CONCLUSIONS 

(a) A feasible computational strategy for the transient solution of dynamically loaded 
bearings or dampers with oil hole feed has been outlined, relying on interpolation 

from previously computed data. ^ n 


<b) 


(c) 


The procedure allows for any number of oil holes, allows for the provision of 
check valves to avoid pump out and allows for oil flow restriction between the 
pply reservoir and the oil hole, albeit at the expense of an added iteration 
The pressure contributions due to the oil holes may be conveniently obtained 

using a cosine power series solution. y ODtamcd 


(f) 


(d) Such solutions are eccentricity and oil hole location (relative to line of centres) 

dependent at any instant of time. ^ ' 

(e) The solutions prove that there is no significant interaction between the pressure 
contributions of adjacent holes, provided they are spaced at least four di hole 
diameters apart. This would pertain in most if not all practical applications. 

The solution procedure is here restricted to the oil holes being located in the 
central circumferential plane. For off-centred oil holes, the series would need to 

difficulties 6 tCrmS 35 Wel1 ’ ThlS eXtenSi0 " iS n0t expecfed t0 create an y un due 
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